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ABSTRACT 

In this work we introduce two new polynomial parametrizations of dark energy and 
explore their correlation properties. The parameters to fit are the equation of state 
values at z = and z — 0.5, which have naturally low correlation and have already 
been shown to improve the popular Chevallier-Polarski-Linder (CPL) parametrization. 
We test our models with low redshift astronomical probes: type la supernovae and 
baryon acoustic oscillations (BAO), in the form of both current and synthetic data. 
Specifically, we present simulations of measurements of the radial and transversal 
BAO scales similar to those expected in a BAO high precision spectroscopic redshift 
survey similar to EUCLID. According to the Bayesian deviance information criterion 
(DIG), which penalizes large errors and correlations, we show that our models perform 
better than the CPL re-parametrization proposed by Wang (in terms of z = and 
z = 0.5). This is due to the combination of a lower correlation and smaller relative 
errors. The same holds for a frequentist perspective: our Figure-of-Merit is larger for 
our parametrizations. 
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1 INTRODUCTION 

More than a decade ago it was discovered that the uni- 
verse is expanding with increasing velocity, and this fact was 
brought to hght by type la Supernovae (SN) observations 
( Perlmutter et al.|199"9l|Riess et al.|1998[[2004}|Astier et al 



|2006 ). Nowadays the accelerated expansion of the universe 
stands as confirmed by several independent observations as 



the already mentioned SNela ( 


Hawkins et al. 


2003 


iGold- 


stein et al.|2003||Readhead et a 


l.|2004||Tegmar 


i et a 


l.|2004| 


Spergel et al.|2007l; the measurements of cluster properties 



as the mass, the correlation function and the evolution with 



redshift of their abundance (Eke et al. 1998 Viana et al 



2002} [Bahcall et al.|[2003l [Bahcall fc Bode||2dd3|); the opt 



cal surveys of large scale structure (Pope et al. 2004 Cole 



et al. 2005 Eisenstein et al. 20051; the anisotropies in the 



Cosmic Microwave Background (CMB) (de Bernardis et al 



2000 Spergel et al. 2003 Komatsu et al. 20111; the cos 



mic shear measured from weak lensing ( Van Waerbeke et al 
2001 Refregier 2003 1 and the Lyman— a forest absorption 
I Croft et al.||1999| [McDonald et al.[2005l ). 

This accelerated expansion can be explained postulat- 
ing the existence of a "new" component in the universe. 
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called "dark energy" which would represent about between 
the 70% — 75% the of the mass-energy of the universe and 
would counteract the effects of the gravitational attraction. 
The remaining 25 — 20% would be composed by dark matter 
mainly and about an extra 5% which would account for the 
baryons, radiation and neutrinos. Although the existence of 
the dark energy is well established observationally, there is 
not a entirely compelling theoretical model which gives a 
theoretical physical framework within which the dark en- 
ergy can be understood, and is of course compliant with all 
main observations. 

Observational data agree rather well with the simplest 
model, ACDM, which puts down the acceleration to the 
presence of a cosmological constant A. This model has be- 
come the standard model of cosmology due to its simplicity 
and compliance with the data, but there are still some points 
that it can not account, e.g. the small value of the cosmo- 
logical constant, which cannot be explained in terms of any 
of the known interactions. 

This situation leads to the proposal of other settings 
which admit a slightly time-variable dark energy as the 
agent producing the the cosmic acceleration. Many such 
models have emerged along theoretical avenues: quintessence 
1999^, Chap lygin gas ([Bento et al.|2002[ ), mod- 



( Zlatev et al. 



ified gravity (Nojiri & Odintsov 2006 



ergy (Li 20041, braneworld models 



holographic dark en- 



Maartens 20051, f{R) 
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theories (Sotiriou & Faraoni 



mensions (Lechtenfeld at al. 



20101, theories with extra di- 



20021, and quite a few others. 



Unfortunately none of them have emerged as definitive an- 
swers for the dark energy problem, so currently the situation 
is one of quick advances on the empirical front (regarding the 
quantity and quality of the data), whereas theory is some- 
what mired in a jungle of alternatives. Since it is not fully 
clear what the observations should be compared with, a try 
and tested approach, borrowed from many areas of physics, 
is defining a parametrization of dark energy functions which, 
if well designed, allows to encapsulate all the observational 
information in a few numbers which can afterwards be com- 
pared with theoretical predictions. 

The most natural approach in this case is considering 
the dark energy equation of state. It is usually assumed that 
this quantity varies slowly with redshift and can be approx- 
imated by a fitting formula with a small number of free pa- 
rameters. These parameters can be constrained comparing 
the ansatz with observations using an optimization proce- 
dure. Choices are typically built upon intuition and prior 
information, but there is yet plenty of room for discussion 
and improvements. This is precisely the direction we follow 
in this work. 

In section [2] we propose two new polynomial 
parametrizations of dark energy. We perform cosmological 
tests to compare them to the most popular dark energy 
parametrization in the literature, both in its classical fashion 
and a new, improved reformulation. We resort to several cri- 
teria to perform our comparison: in addition to the value, 
we also pay attention to the correlation of the coefficients, 
figure of merit, fractional errors on the free dark energy pa- 
rameters and the deviance information criterion (DIG). The 
reason for an analysis from several perspectives is that we 
find to be an insufficiently informative tool, as it does 
not offer any reward upon some important improvements 
like tighter constraints or lower correlation. 

Our tests make use of Type la supernovae (SN) and 
baryon acoustic oscillations (BAO), which are so far the best 
representatives in the categories of standard candles (ob- 
jects with well determined intrinsic luminosity) and stan- 
dard rulers (objects with well determined comoving size). 
Such probes provide us with distance measures related to the 
Hubble factor H{z), and as they are low redshift datasets, 
their suitability to constrain dark energy is strong, as this 
component of the cosmic soup has has began to govern the 
evolution of the Universe just recently, according to most 
evidences. This combination has additional advantages: SN 
measurements come in the form of luminosity distances and 
therefore give a smeared out information on H{z) (two in- 
tegrations are required) but still remain extremely useful 
because of the large number of measurements and their con- 
siderably good quality. On the other hand, BAO measure- 
ments, though currently scarce, involve 1/H{z) directly, so 
they are expected to favour sensitivity considerably, besides 
being of even better quality than SN data. In principle, one 
could also consider including CMB data, but typically they 
do not improve constraints on dark energy parameters signif- 
icantly (WMAP7-year data alone constrain w in quiessence 
models with about a 40% error), and the dark matter den- 
sity Q,rn is generally the only parameter on which those data 
exert a strong impact. A good compromise between simplic- 



ity and advantages offered by CMB as regards Sim is the use 
of priors, and that is the approach we use. 

Interestingly, we will not only use the latest observa- 
tional data to obtain the constraints on the parameters of 
the models, but we will also consider mock data simulating 
a forthcoming survey, as described in sections |3] We present 
the mock datasets for the two main measurable quantities 
expected from a line-of-sight, high-resolution spectroscopic 
baryon acoustic oscillations survey ( Percival|2010 EUCLID 
EUCLID! [Refregier et al.|p010l [Beaulieu et al.||2010|, asl t 
has done in a previous work ( Escamilla- Rivera et al. 112011 1. 



We combine these data with synthetic pre-WFIRST 
(Wide-Field Infrared Survey Telescope) supernovae data to 
throw further light on the constraining power and suitabil- 
ity of the parametrizations proposed, so we are allowed to 
strengthen our conclusions. 



2 DARK ENERGY PARAMETRIZATIONS 

The Friedman equations explain how an homogeneous and 
isotropic universe expands in the context of General Relativ- 
ity (or generalizations). Considering such a universe is filled 
with several fiuid components with pressures pi and energy 
densities pi, those equations will read 



8nG 



k 



inG 



^(p, + 3pi 



(1) 



(2) 



where a and H are the scale and Hubble factor respectively. 
At present, the two main such fluids are dark matter and 
dark energy, the last being the governor. Dark energy (de) 
enters our picture in an effective way in the sense we do not 
appeal to any fundamental theory to deduce its behaviour 
and just let it be represented by a phenomenological equa- 
tion of state (EOS): 



w[z) — 



pde 



(3) 



The importance of the EOS is significant, because it de- 
termines the form of the Hubble parameter H{z) or any 
derivation of it which is necessary to obtain the observable 
quantities. Under the spatial flatness assumption (fe = 0) 



= n^{l + zf + naeX{z) 

with 



X{z) 



Pdejz) 
Pde(O) 



exp 3 



1 + w{z) 

1 + 2 



dz 



(4) 



(5) 



and Qde = 1 — ^m- 

Considering the huge number of contributions to the 
topic, and the knowledge so far gathered it might seem hard 
to make improvements, but we believe some avenues opened 
by Wang in (Wang 2008^ are worth exploring and allow for 
little but valuable advances towards parametrizations that 
make the best out of the data available. Questions like the 
most convenient choice of the two dark energy parameters 
or the most suitable model comparison criterion are worth 
being looked at once and again, particularly to build a more 
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solid background to make the most out of future data and 
their expected far better statistical value. 



2.1 Chevallier-Polarski-Linder parametrization 

The Chevallier-Polarski-Lindcr (CPL) parametrization, first 
discussed i n (jChevallier fc Polarski_ 200T] ) and reintroduced 
in ( Linder|[2003 1, defines the dark energy equation of state 



w{z) =Wo+Wa 



1 + 2' 



(6) 



where wo is the value of the dark energy equation of state 
today (i.e. at redshift z — 0). This parametrization has been 
widely used because of its simplicity, sensitivity to observa- 
tional data and because it is well behaved and bounded at 
high redshifts. Its adoption by the Dark Energy Task Force 



I Albrecht et al. 2009 1 as a preferred parametrization has 



contributed to its popularity. 

However, this parametrization has shortcomings: its 
fiexibility to determine cosmological parameters for some 
dark energy models with rapid evolution has been put in 
doubt recently, and on the other hand the second parame- 
ter is typically very poorly constrained thus losing power of 
conviction about the conclusions to be drawn from it. 



2.2 Wang parametrization 

The quest to delineate the expansion history of Universe is 
expected to make big ambitious moves in the future. As a 
result it is expected that many more and better observa- 
tional data will be available. Obviously it is worth getting 
prepared to making the best profit out of this avalanche of 
data to come, and this can be done by learning as much as 
possible from the data we already. Intuitively, in order to be 
able to discern which dark energy model is closest to reality 
and to offer a sensible interpretation of the results it seems 
necessary to consider dark energy parameters which offer 
clear advantages. Apart from having a clear physical mean- 
ing, for a parameter to be eligible there should be previous 
hints or theoretical grounds forecasting reasonably tight con- 
strains, and as it is well know, the parameter should be the 
least possibly correlated to others. Suitability criteria for 
parameters along those grounds can be sketched with the 
help of comparison tools rewarding for those nice features. 
The (frequentist) FoM favours low correlation whereas the 
(Bayesian) does not only favour that feature, but it also re- 
wards for tight constraints also, so both these criteria (FoM 
and DIG) do poorly when high correlation is present in the 
dark energy parametrization one considers. In this respect, 
one must make a fair use of those statistical tools. For in- 
stance, if one takes performs a reparametrization of a certain 
scenario, the FoM and DIG as calculated for both cases will 
typically be direct indicators of improvements (or worsen- 
ings) in terms of correlation, but changes in the value of 
those quantities should not be held in the same grounds 
as the changes occurring when constraining the same single 
parametrization with two different datasets. 

The GPL parametrization suffers from quite a signifi- 
cant correlation, besides the fact that constraints on one of 
its parameters are typically large in percentual terms. But 
on the other hand it has some nice features, so a convenient 



redefinition as concerns those two issues was put forward to 
yield an improved situation, although the encoded informa- 
tion remains exactly the same. As mentioned, this was done 
m fWang..2008j , where a new dark energy description was 
given in terms of its value at present, wq, and at redshift 
z = 0.5, ii;o.5, was given. Explicitly 



w{z) — 3iuo.5 — 2it)o + 



3(uio — TOO. 5) 
l^z ' 



(7) 



was proposed, and this can be seen just as a rearrangement 
of the classic GPL parametrization. This reformulation mi- 
norates the correlation between the parameters and allows 
to obtain tighter constrains along with a more transparent 
interpretation of the parameter estimation results, just in 
the spirit of making an optimal use of observational data 
from future surveys. 

The results in our paper reinforce the view that Eq. [7] 
is a preferred way of exploiting the GPL parametrization, 
and on the other hand we suggest, with the help of two 
new parametrizations, that the (wcwcs) couple is indeed 
a very good choice despite the specifics of the dark energy 
evolution, provided it is smooth enough and with bounded 
early asymptotic behaviour. 



2.3 Polynomial parametrizations 

The two proposals we are making are somewhat inspired on 
the one hand by the GPL parametrization and on the other 
hand by an a proposal consisting in a expansion in powers 
of the quantity {1 + z) which emerged naturally from the 
relationship between the redshift and the scale factor and 
showed computationally convenience. This second inspiring 
setup we are referring to was proposed for the first time in 
( [Weller fc Albrecht|2002l ) in the form 



w(z) 



-1 + Ci{l + Z) + C2{1 + Z) 



(8) 



However, this sort of parametrization poses problems at high 
redshifts, as |to(2:)| grows unboundedly with z and one will 
either end up with a superphantom model or a superluminal 
one. This motivates considering generalizations which are 
devoid of this pathology and may match or event surpass 
the nicety of GPL or its reformulation (recall Eq. [7|. 

Two possible routes that retain some similarity with the 
previous case but also some improvements arise from 



W{Z) = -1 + Ci(l + f{z)) + C2(l + f{z)f 



(9) 



with f{z) a smooth and at the same time simple function or 
the slightly more general 



W{z) = -1 + Ci ffi(l + f{z)) + C2 92(1 + .f{z)), 



(10) 



where gi and g2 are some smooth and simple functions as 
well. 



2.3.1 Conventional polynomial 

In the spirit of the first scheme above we propose 

= (11) 
1 + z 

so we avoid high-redshift unboundedness by the same via as 
in GPL. Therefore we have 
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w{z) = -1 + Ci 1 + 



1 + z 



C2 1 + 



1 + Z 



(12) 



which we dub conventional polynomial parametrization. If 
more compactness is desired one can also write the latter as 



w{z) = — 1 + Ci 



1 + 2z 
1 + z 



C2 



1 + 22 
1 + Z 



(13) 



It is convenient to leave out the constant term of value —1 
from the computational point of view because as one does 
not expect a large departure from a ACDM setting, it is 
reasonable to confine (at least initially) the parameter search 
region to |ci| < 1 and to \c2\ < 1. 

However, as we have discussed already it is desirable 
to fit parameters which are more or less physically trans- 
parent and it in addition are just lightly correlated. At a 
first stage we wish to compare our conventional polyno- 
mial parametrization with Wang's, so the best way to do 
so is to consider exactly the same two dark energy param- 
eters. Then, at a later stage, we will check whether the 
good behaviour as correlation is concerned is shared by our 
parametrization. With those arguments in mind we refor- 
mulate the proposal made in Eq. [l3] by letting 



ci = -(16it)o - 9uio.5 -I- 7), 



C2 — —3wo + 



(14) 
(15) 



This way, we move on to a scenario in which wo and wo.s are 
the parameters subject to estimation, which can be explicitly 
reconstructed using Eq. 



X(Z) = (l + ;2)f(-8'"0+9»0.5 + l) 



(16) 



which consistently comes down to the ACDM case for wo ~ 
Wo. 5 = — 1- 

For completeness and purposes related to model selec- 
tion, it is convenient to compute and effective Wa parameter 
for this model. By analogy with the CPL case we define 



Wa = lim w(z) — Wo, 

z—^oo 

so in this case 

9 1 

Wa = — OUIQ + 2'^0.5 — 2 

is the result we get. 



(17) 



(18) 



2.3.2 Chebychev polynomial parametrization 

Now we want to make a further generalisation in the fash- 
ion of our general proposal above, by considering a bit more 
involved functions. In this case we make use of Chebychev 
polynomials (of the first kind), which have a significant role 
in most areas of numerical analysis, as well as in other areas 
of Mathematics (polynomial approximation, numerical inte- 
gration, and pseudo spectral methods for partial differential 
equation, etc.) 

Once again we have to take into account that we can not 
constrain accurately more than two dark energy parameters 
( [Linder &: Huterer|2005[ ), thus we have cut the expansion at 
the second order. Specifically we propose 



W{z) = -1 + CiTi(l + f{z)) + C2T2{1 + f{z)), 



(19) 



with Tn being the first kind Chebyshev polynomial of degree 
n and f{z) — z/{l + z) as before. A convenient presentation 
of that parametrization is 



^ ^ 1 , / 1 + 2z 
w(z) = -1 -I- ci ( ' ^ ^'^ 



l + 2z 
1 + z 



(20) 



In this case too we switch to more amenable parameters and 
then let 



ci = —(23^0-9^0.5 + 14) 
3 

C2 = -^^(4^0 - 3uio.5 + 1). 



(21) 
(22) 



Then, one should simply resort to Eq. ([5| to produce the 
whole scenario: 



X{z) = (l + 2)-i 



|(52u,o-45™o.5+7) 



3z(rao('>S^ + 56)-9™o 5 (6^ + 5) + 14z+ 1 1) 



(23) 



For this first attempt at depicting a dark energy dom- 
inated universe with Chebyshev polynomials in terms of 
those functions we choose wo and wo.5 as our parameters, 
but it it would be not surprising than one could do better if 
some extra work was done in the direction of making param- 
eter correlation smaller. For the time being we just pursue to 
compare directly our polynomial proposal with a preferred 
presentation of the CPL parametrization on the one hand 
and with our conventional polynomial proposal on the other 
hand. 

In that direction and as it has been done before, we 
compute an effective Wa parameter which in this case takes 
the form 

Wa = ^ (-49W0 + 45u;o.5 - 4). (24) 



2.3.3 Some additional remarks 

In the first place, these two new routes can be viewed as 
perturbations of ACDM, particularly, at low redshifts, which 
is the region most accurately described by the current data. 
If one wants to perturb ACDM in either the CPL or the 
Wang scenarios wo has to be anchored at —1 and then there 
is only one free parameter to play with. In contrast, our 
two new models may model two-parametric departures from 
ACDM, and thus have more flexibility in principle. 

Note that it would be possible to consider the two new 
parametrizations along with CPL(Wang) if one let the pa- 
rameter space have one more dimension. Indeed, if we let 
our parametrization be of the form 



w{z) = 61+62 



l + 2z 
1 + z 



+ bs 



2z 



1 + z 



(25) 



then our conventional polynomial case would be obtained 
for hi — —1, 62 = cl, and 63 — C2; the Chebyshev case 
would follow from the choice 62 = ci, 63 = 2c2 and 61 = 
— (1 + C2); and finally the CPL case would be obtained from 
61 = wo — Wa, 62 = Wa and 63 = 0. 
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3 OBSERVATIONAL DATA 

As dark energy is expected to have started to dominate at 
recent times, low redshift datasets are the obvious choice to 
put the tighter constraints on each dynamics, whereas high 
redshift ones may be viewed as complementary. Thus the 
combination of SN and BAG datasets, given their quality in 
both cases, and the considerable number of data points in 
the case of the SN, is an excellent choice given the state of 
the art. Besides, new avenues on BAG ^Bassett & Hlozek 
[20I0I ) are to be open soon which will allow to exploit the 
tremendous potential of this new astronomical tool towards 
constraining the main evolutionary features of dark energy. 
As we have already mentioned, one of our objectives is to 
introduce new promising parametrizations as alternatives to 
one of the commonest, but one of the other objectives is to 
forecast how the old parametrizations and our challengers 
will cope with new data. 

The literature provides a large number of papers where 
simulated supernovae data are used in the way we have just 
mentioned, but to our knowledge synthetic baryon acoustic 
oscillations data have only been presented and exploited in 
I Escamilla- Rivera et al. 20111. This builds on considerable 
theoretical efforts in different forecast aspects (Blake ct al.l 



2006), which have crystallized in the package Initiative for 
Cosmology (iCosmo) (see (Refregier et al. 2011 1) and its 



BAG modules, which have allowed us to produce these mock 
data. We have modified and extended this general purpose 
software to produce mock data from a line of sight, high- 
precision BAG spectroscopic survey as the one described in 
( Percival|20l'0l|EUCLID|EUCLID[ ) and pre-WFIRST super- 
novae data. 



3.1 Baryon Acoustic Oscillations 

Baryon Acoustic Gscillations (BAG) have emerged as a 
promising standard ruler in cosmology, enabling precise 
measurements of the dark energy parameters with a min- 
imum of systematic errors ( |Wang et alT][20Tol |Blake et al.] 
20061. These oscillations were originated before recombi- 



nation due to the density fluctuations created by acoustic 
waves generated by primordial perturbations. After recom- 
bination, photons decoupled and propagated freely leaving a 
signature of the primordial perturbations in the CMB tem- 
perature distribution. A similar but attenuated feature ap- 
pears in the clustering of matter; the peaks and troughs of 
the acoustic waves gave rise to overdense regions of baryonic 
matter which imprint a correlation between matter densities 
at the scale of the sound horizon at recombination 



,{Zr) = / 



H(Z) 



dz, 



(26) 



where Cs is the sound speed (Bassett & Hlozek 2010 Per 



cival 20061. Cosmological Microwave Background (CMB) 



anisotropies provide an absolute physical scale for these 
baryonic peaks, but these features can be also determined 
in a galaxy survey. By comparison of the absolute value 
given by CMB and the observed location of the peaks of the 
two-point correlation function of the matter distribution, 
^{z) given by a galaxy survey, one can obtain measurements 
of cosmological distance scales. Measuring the BAG scale 
from galaxy clustering in the transverse and radial directions 



yields measurements of r{z) /rs{zr) and of rs{zr)H{z) (Blake 



fc Glazebrook|[2003l |Hu fc Haim"anl[2003l [Seo fc Eisenstein 
20031), where 



cdz' 

/o 

is the comoving distance at redshift z. 



3.1.1 Percival et al. 



(27) 



In ( Percival W. et al.|2010 ) Gaussian values on the distance 
ratio, r3{zdia.g)/Dv{z), at redshifts z — 0.2 and z = 0.35, are 
given from the measures obtained by combining the spec- 
troscopic Sloan Digital Sky Survey (SDSS) and the Two- 
Degree Field Galaxy Redshift Survey (2DFGRS) data. This 
distance ratio represents the comoving sound horizon at the 
baryon dragging epoch, Zdrag, 



(^drag ) 



H{z) 



dz 



(28) 



over the effective distance Dv{z), defined as (Eisenstein 
et al.|2005| ) as 



Dv{z) = 



(l + zfD\{z) 



H{z) 



1/3 



(29) 



Da being the angular diameter distance which takes the 
form 



Da{z) = 



1 + z 



Az' 
Hiz'Y 



(30) 



in a flat universe containing only matter and dark energy. 

In order to estimate the dark energy parameters in the 
context of Bayesian statistics, as detailed in Appendix |Aj 
we need a definition of the which reflects the difference 
between the observational data and the values given by our 
models. In our case this requires giving an expression for the 
comoving sound horizon at the baryon dragging epoch, and 



we have used the fitting formula proposed in (Eisenstein & 
Hu|1998| ): 



rs(2drag) = 153.5 



0.02273 



0.1326 



(31) 



Now, taking into account the Gaussian values at 2; = 0.2 



and 0.35 from the BAG data in ( jPercival W. et al.|[2010| |, 

we can calculate 



Xbao = 
where 

fs (^drag , 



^ BAOn / ^ 



-l-iBAO/ 



flrn,i^b', 0) rs{zdrag,^m.,i^b', 0) 



Dv{0.2,Q,n;e) ' L>v(0.35,a^;6>) 



^BAO ^ (0.1905,0.1097) 



and 



30124 



17227 



-17227 86977 
being the inverse of the covariance matrix. 



(32) 



(33) 



(34) 



(35) 
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3.1.2 High precision spectroscopic redshift surveys 
Future BAO surveys are expected to represent a real break- 



through in our knowledge of the dark energy. In (Percival 
2010 EUCLID||EUCLID I is proposed a survey to measure 
spectroscopic redshifts for 6.1 x 10^ luminous galaxies and 
clusters of galaxies out to redshift z — 2 over 20000 deg^ , 
reaching a dz/{l + z) < 0.001, enough to resolve the BAO 
feature along the line of sight, and achieving much better 
dark energy constraints than their predecessors. 

We have used the Initiative for Cosmology (iCosmo) 
software package to generate the BAO mock data. In this 
case are concerned with the following two signatures of the 
BAO peak: 



y{z) ■ 



y (z) = 



rs{Zr) ' 

r'{z) _ 



(36) 



(37) 



_ ^/H{z)^ 
raiZr) rs{zr) 

The publicly available code has built-in routines based on 
the universal BAO fitting formulae for the diagonal errors 
on y and y' presented in ( Blake et al.||200^ K Once we have 
made the proper modifications of the code to replace the de- 
faults (Peacock and Dodds (Peacock & Dodds|l996) power 
spectrum and Small et al. ( jSmail et al.|1994 1 galaxy distri- 
bution) with the survey properties (those quoted above and 
besides Zmed ~ 0.46, redshift range 0.1 < z < 0.9), we have 
written extra codes to generate a large number of normal 
random realizations around a fiducial model which is specif- 



ically the wcdm-|-sz-|-lens case from WMAP7-year ( [Komatsu 
|et al.|201l| , which has n„ = 0.26 ± 0.099 and is phantom- 
like with w = —1.12 ± 0.43., Then after some reduction the 
synthetic BAO data presented in tablelllhave been obtained. 
In addition, in the the corresponding we have introduced 
priors based on the values of the matter and baryon density 
presented in ( Burigana et al.||2010 l as a forecast analysis of 
Planck: using the result Qrnh'^ = 0.1308 ± 0.0008 we con- 
struct a weak Gaussian prior, whereas with Qbh^ = 0.0223 
we construct a fixed prior. In both cases we use h — 0.742 as 
is given by ( [Riess et al.|2009[ ). See as well Fig. [l] for a graph- 
ical account of the features of our BAO simulated data. 

In this case we need an expression for the sound horizon 
at recombination to constrain the dark energy parameters: 

r,(a,) = I , M/2 M , oM/2 ■ (38) 



n 



(a + ae,)i/2 (i + i?)i/2 



which can be evaluated as described in (Efstathiou & Bond 



1999 1 and gives 



rs{ar) - 



4000 



(39) 



hi 



[1 + R{z 



11/2 



+ [R{z,) + R, 



11/2 



Mpc 



where r^i, = 0.6813 is the ratio of the energy density in neu- 
trinos to the energy in photons. The parameter R = 3pb/4:pry 
is numerically obtained 

R{a) = 30496!r26/i^a. (40) 

The scale factor at which radiation and matter have equal 
densities is 

1.6813, 



a^g' = 24185 



1 + 



(41) 



and the redshift of recombination Zr is given by (Hu & 
Sugiyama|1996 1 in the following fitting formula 

Zr = 1048 (1 + 0.00124(!r26/i^)~°-^='*) (1 -I- gi (fi^/i^)"^), (42) 

being 



(43) 



(44) 



gi = 0.0783{nbh'^y°-^^^ (1 + 39.5 {nth'^f-''^^ 
and 

52 = 0.560 (1 + 21.1(f7i,/i^)^ " . 

The function for BAO mock data is now defined in 
Eq.|45], where Nmode is the number of mock data, in this case 
4. And 6 = 61,62, ... are the dark energy model parameters. 
Here we have conveniently accounted for the slight degree 
of correlation existing between y and y' , and as suggested 
in ( |Seo fc Eisenstei"n||2007[ ) we will fix for our calculations 
Py,y' = 0.4. 



3.2 Type la Supernovae 

Type la supernovae (SN) are the explosions that take place 
at late stages of the stellar evolution. They have been recog- 
nized as a powerful probe of cosmological dynamics, as they 
give a good measure of the cosmological expansion rate. Su- 
pernovae provided the first probe for the cosmological ex- 
pansion ( Riess et al.||1998 Perlmutter et alT]|1999p and are 
considered as standard candles ( Leibundgut|200T \ . 

The statistical analysis of such SN samples rests on the 
definition of the modulus distance: 



p{zj) = 5log^a[d'L{zj,e)]+ fio, 
where dL{zj,6) is the luminosity distance: 

1 



dL{z,e) = (1 + ^) 



dz'- 



H{z',0) 



(46) 



(47) 



The best fits for the parameters of a given model can be 
obtained by minimizing 



j=i 



{fith{zj; fio, 0) - ^.ohs{zj)Y 



(48) 



glected correlations between measurements at different red- 
shifts, as they are typically small; as it is well known, doing 
this just induces a slight degradation in the constraints, and 
therefore it is simply a conservative and acceptable proce- 
dure. 

In our xsN we have a nuisance parameter, /io, which 
encodes the Hubble parameter and the absolute magni- 
tude M, and has to be marginalized over. However, when 



one works with an homogeneous data sample (Nesseris & 



Perivolaropoulos 20041, an alternative method is used for 



that purpose. This method consists in max imizing the like- 
lihood by min imizing with respect to /io ( Elgaroy & Mul- 
tamaki|2006 |. Then one can rewrite the Xsn(a*Oi 0) as 

(49) 



Xs-n(6) = ci - 2c2^o + C3M0 



being 



Cl 



3=1 



(p.ohs{zj) - 5 logio dL (zj;0)) 



(50) 
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XbaoW ■ 



V'V L j=i 



y ,3 



2 Py,y' ^ ^ 



' {y{Zj,VL^,9.b\0) - ymock{Zj))iy'{zj,nmMb;d) -y'mockiZj)) 



^yj'^y',] 



(45) 



Table 1. Mock BAO data from 



0.6 15.005 ±0.117 21.401 ±0.291 



0.8 
1.0 
1.2 
1.4 
1.6 
1.8 
2.0 



19.221 ± 0.109 
22.708 ±0.102 
25.762 ±0.098 
28.635 ±0.096 
31.126 ±0.101 
33.540 ±0.107 
35.357 ± 0.112 



18.953 ±0.186 
17.012 ± 0.137 
15.098 ±0.101 
13.421 ± 0.079 
12.213 ±0.069 
10.987 ± 0.061 
10.094 ± 0.055 



-s-^ Mobs (zj) - SlogiprfL {zj;9) 



C3 



(51) 
(52) 



population at z > 1.6. Specifically, the redshift distribution 



With the minimization over /io one obtains an expression in 
terms of C2 and C3: 



Mo = C2/C3, 

so the function takes the form 



xiN(0) = ci-^. 

C3 



3.2.1 Union2 



(53) 



(54) 



Union2 ( Amanullah et al.|2010 | is one of the largest type la 
Supernovae samples up to date, and it consists of 557 SNIa 
and covers a redshift range from to 1.4. This sample has 
increased the number of well-measured Type la supernovae 
of the previous Union ( [Kowalski et al.|2008[ ) at high redshifts 
by the combination of different data sets. The Union Sample 
has been extended with six type la supernovae presented in 
([Amanullah et al.||2010| ), the SNe from ( [Amanullah et al 

2008) , the low-z and intermediate-z data from (Hicken et al 

2009 ) and ( [Holtzman et al.|2008| respectively 



This sample has been obtained after some improve- 
ments in the Union analysis chain: all light curves have been 
fitted using a single light curve fitter (SALT2) to eliminate 
differences and systematic errors have been directly com- 
puted using the effect they have on the distance modulus. 



3.2.2 Mock data for a pre- WFIRST stage 

To create SN mock samples we have reproduced the pre- 



WFIRST observational simulation as reported in (Albrecht 
[et al.|2009[ ), which has two population peaks, one at z < 0.1 
and the other at 0.6 < 2 < 0.7; along with a very scarce 



suggested in ( Albrecht et al. 2009 ) is reproduced here in 
Table H 

The formulae for errors on SN magnitudes that we use 
follows from a prescription used in the binning approach 
( |Kim et al.|2004| [Upadhye et al.|2005| [Ishak|2007[ ) , in which 
they are calculated as follows: 



where: 



1 ^pec 1 ^ syst 5 



(55) 



• o-int = 0.15 is the intrinsic dispersion in magnitude per 
SN, assumed to be constant and independent of redshift for 
all well-designed surveys; 

• o-pec ~ 5o"t,/(ln(10)c2) is the error due to the uncer- 
tainty in the SN peculiar velocity, with = 500 km/s, c is 
the velocity of light and z is the redshift for any SN; 

• o-syst = 0.02{z/zmax) is the floor uncertainty related 
to all the irreducible systematic errors with cannot be re- 
duced statistically by increasing the number of observations. 
The value 0.02 is conservative from the perspective of what 
space-based missions could achieve. Those are precisely the 
resources expected to provide high-redshift SN, which are in 
turn the ones in which the systematic error term is expected 
to contribute. Note as well that Zmax is the maximum ob- 
servable redshift in the considered mission and this linear 
term in redshift is used to account for the dependence with 
redshift of many of the possible systematic error sources (for 
example the Malmquist bias or gravitational lensing effects). 

We have included some extra (though very slight) noise, and 
then we have checked that our mock data are compliant with 
the main features of the Union2 sample. The fiducial model 
used is again the wcdm+sz+lens case from WMAP7-year as 
quoted above , which has Q,m = 0.26±0.099 and is phantom- 
like with w = -1.12±0.43. The function for the SN mock 



data has been constructed as described in Sec. (3.2 1. Then 



we have let iCosmo generate the dh values for the redshifts 
in the table a large number of universes drawn randomly and 
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redshift, z redshift, z 



Figure 1. Fractional errors on synthetic BAO data for y (a) and y' (b). 
Table 2. Redshift distribution of pre-WFIRST SN samples 



redshift bin # of SN 



< 0.1 


500 


0.1 - 0.2 


200 


0.2-0.3 


320 


0.3 - 0.4 


445 


0.4 - 0.5 


580 


0.5 - 0.6 


660 


0.6 - 0.7 


700 


0.7 - 0.8 


670 


0.8 - 0.9 


110 


0.9 - 1.0 


80 


1.0 - 1.1 


25 


1.1 - 1.2 


16 


1.2 - 1.3 


16 


1.3 - 1.4 


4 


1.4 - 1.5 


4 


1.5 - 1.6 


4 


> 1.6 


4 



normally distributed around the fiducial one, and finally we 
have performed a reduction to give our mock sample. 



4 RESULTS 

Following Bayesian Statistics, we have inferred the values 
Qrn and the dark energy parameters for the models consid- 
ered. For further details see |A] We report our findings in 
two main ways: on the one hand we present our best fits, 
errors and derived quantities in Tables [3] [4] on the other 
hand we present credible contours obtained after a numer- 
ical marginalization over flm, where we have approximate 
the likelihood as a Gaussian and then analytically marginal- 
ize, see ( Taylor Sz Kitching|2010 1 and references therein. As 
mentioned, we have considered the combination of real and 
mock SN and BAO data as discussed above, and in addition 
we have introduced in all cases a Gaussian prior on Qrn and 
Qb deduced in (Burigana et al 
sis of Planck 



2010) 



Qrnh 



as a forecast analy- 
0.1308 ± 0.0008, Qt = 0.0223 with 



h = 0.742 as is given by ( |Riess et al.|2009 l. 

The presence of the prior leads to uniformity in the best 
fit value of Qm thus minorating its influence in the dark 
energy constraints. This trick, together with the marginal- 
ization, allow to focus the discussion on the dark energy 
parameters thus offering a clearer picture of what each 
parametrization can offer. 



Let us now examine our results using different criteria. 
First of all uncertainties in percentage terms on the dark en- 
ergy values can be considered. The conventional polynomial 
turns out to give the lowest percentage error on wq for both 
current and recent data, and the second best is the Cheby- 
shev setting, although we only get a marginal difference with 
respect to the conventional polynomial and Wang's model. 
The situation gets reversed, though, for the percentage er- 
ror on wo.5 , and Wang's parametrization is the best per- 
former, but again differences are small when compared to 
our parametrizations (CPL gets excluded from this discus- 
sion item as the parameter wo.s is not considered). 

Note that the substantial reduction on percentual er- 
rors when one moves from current to mock data, typically 
they become almost three times smaller. Up to some degree 
this may be a consequence of the use of a fiducial model in 
the simulated data, probably real future data will be signif- 
icantly better than current ones but perhaps not as remark- 
ably as our synthetic data. 

Summarizing, this section of the analysis, we see that 
in general choosing wo and wo.s as the parameters to con- 
strain is worthy as percentual errors are low. In our two 
new parametrizations are indeed valuable from this particu- 
lar perspective, but are they as valuable when one considers 
other criteria? 

Another interesting point of view to interpret our re- 
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0.5 1.0 
redshift, z 

Figure 2. Mock data obtained with iCosmo for the pre-WFIRST specifications. 



Table 3. Constraints on darlc energy parameters and derived quantities from current data. 



Model 






Dark energy parameters FoM^ang Pi, 2 DIG 








(FoMdetp) 



CPL 


544.91 


0.315 ±0.033 


wo = 


-1.033 ±0.180 


14.16 








Wa = 


-0.742 ± 1.520 


(14.16) 


Wang 


544.91 


0.315 ±0.033 


Wo = 


-1.033 ±0.180 


42.47 








WO. 5 = 


-1.281 ± 0.361 


(14.16) 


Chebychev Polynomial 


544.96 


0.314 ± 0.032 


WO = 


-1.049 ±0.163 


43.63 








«'0.5 = 


-1.268 ±0.373 


(11.67) 


Conventional Polynomial 


544.98 


0.314 ± 0.032 


Wo = 


-1.055 ±0.157 


44.11 








WO. 5 = 


-1.263 ±0.377 


(9.80) 



-0.930 10.30 



-0.792 8.11 



-0.785 8.02 



-0.780 7.98 



suits is Pearson's correlation coefficient (for the dark en- 
ergy), defined as 



As explained in ( Wang 2008 1 , this specific definition has 



Pi, 2 ~ 



„2 
^^12 

0"lO"2 



(56) 



which can be used to study the lineal correlation between 
either two dark energy parameters. Here ai2 stands generi- 
cally for the non-diagonal element of the covariance matrix 
for the parameters 1 and 2. A value of pi,2 close to zero 
will tell us there is no linear correlation between them, and 
lowering correlation typically improves constraints, so if a 
given parametrization achieves that goal naturally then it 
will most likely provide an overall worthy scenario. 

A related magnitude is the (frequentist) Figure-of-Merit 
(FoM), which has been defined in two slightly different but 
related ways in the literature. We present both and discuss 
their specific features and report their values for the benefit 
of the readers keener to one version or the other. Following 
the corresponding references we use subscripts to differenti- 
ate them. According to ( Wang|2008 1 

FoMwang = ^====^=^=, (57) 

^/det C(ci, C2, C3, ...) 

where C(ci, C2, ca, ...) is the covariance matrix of the corre- 
sponding Ci dark energy parameters one is concerned with. 



two advantages: firstly it is easy to calculate for either real 
or simulated data, and secondly it has an easy physical in- 
terpretation. The FoM penalizes experiments that yield very 
correlated estimates for the dark energy parameters. Hence 
the FoM is larger when the dark energy parameters Ci are 
chosen such that they are minimally correlated with each 
other. But obviously, if the model considered does not have 
a low degree of correlation per se one will be giving a poor 
estimation of the performance of the observational tests. 

A closely related, and perhaps more popular version of 
the FoM was proposed by the "Dark Energy Task Force" 
(see ( [Albrecht et al.|2009[ ) for a recent revision of the topic), 
and it simply reads 



FoMe 



y/det C{wo, 



(58) 



Thus, the main difference with respect to the FoMwang is 
that in this view the two relevant parameters are wo and Wa, 
and not any others. However, Wa is a parameter informing 
about early asymptotics, which is a region with extremely 
hard observational access and physical interpretation, and 
thus, the debate persists of whether a FoM using that pa- 
rameter may not be given artificially large results. 
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Table 4. Constraints on dark energy parameters and derived quantities from simulated data. 



Model 






Dark 


energy parameters 


FoMwang 

(FoMdetf) 


Pl,2 


DIG 


CPL 


5320.38 


0.269 ±0.005 


Wo 


= -1.151 ± 0.041 
= 0.244 ± 0.207 


616.25 
(616.25) 


-0.946 


21.44 


Wang 


5320.38 


0.269 ±0.005 


Wo 


= -1.151 ± 0.041 
= -1.069 ±0.031 


1848.70 
(616.25) 


-0.712 


11.69 


Cliebychcv polynomial 


5320.43 


0.269 ±0.005 


WQ 


= -1.140 ±0.035 
= -1.073 ±0.029 


2065.28 
(504.85) 


-0.666 


11.38 


Conventional polynomial 


5320.45 


0.269 ±0.005 


WQ 


= -1.137 ±0.033 


2147.57 


-0.648 


11.29 



«,o.5 = -1.075 ± 0.028 (477.24) 



According to Appendix [C] the two FoMs are linearly 
related: 



FoMdetf = msFoMwang 



(59) 



with 7713 = 1/3, 2/9, 11/45 for Wang, Chebychev polynomial, 
Conventional polynomial models respectively. 

Our results on the FoM get summarized very simply: 
the conventional polynomial has the largest value of the 
FoMwang, whereas the results get complete reversed for the 
FoMdetf - In addition, the ratios between the different FoMs 
are very similar for real and synthetic data. However, as fol- 
lows from our discussion before, the price paid to get a larger 
FoM (by a redefinition) it to waive the importance of cor- 
relation, and further investigations and reflections would be 
needed to provide a definite solution to this debate, which 
is, on the other hand, out of the scope of this paper. 

For both sorts of data the conventional polynomial 
parametrization is naturally less correlated than all three 
others. The second best is the Chebyshev one, and all three 
are considerably less correlated than CPL. Nevertheless, 
even though Wang's second parameter Wc = wo.5 is chosen 
for the low degree of correlation with wq , a better choice for 
that purpose would be the dark energy EOS parameter eval- 
uated at a lower redshift (say z ~ 0.25). This was already 
pointed out i n (|Wang||2008| ) and confirmed in ( |Escamilla-| 
Rivera et al. 20111. In agreement with our discussion the 



degree of correlation at a certain low redshift in our two 
new parametrizations immediately drives to very narrow er- 
rors on the total w at that location, as reflected in Fig. |3j 
where we represent the total w for our best fits and the 
Icr error bands. In fact these two figures serve the addi- 
tional purpose of illustrating the sort of evolution described 
our parametrizations. Nevertheless, as correlation is a topic 
worth of further consideration we elaborate in Appendix [B] 
Finally, the last criterion we resort to is the Bayesian 
deviance information criterion (DIC). This is a very inter- 
esting way to examine results in this context, in particular 
when one finds marginal differences in the values. This 
criterion accounts for the dependence of on our param- 
eters, but it also somehow involves the correlation and the 
percentual errors, thus turning out to be far more informa- 
tive and having more discerning power. When applied for 
model selection, as we put forward, the setting with the 



lowest DIC is in principle the best. In our case, as we are 
mainly concerned with dark energy issues given that Q.m is 
very tightly constrained we construct our DIC starting from 
a marginalized over Q,m- Namely 



DIC = 2x^(0) -X^(06/) 
where 



-2 log ( j exp(-x72)dn„ 



(60) 



(61) 



The behaviour of the DIC follows the same pattern as 
the FoM (although the FoM is not so convenient as a model 
selection criterion as it does not involve y^). Basically, the 
conventional polynomial model is the best one, then we have 
the Chebyshev polynomial model, then Wang's scenario, and 
finally, the CPL model closes the ranking with the highest 
DIC by far. 

Finally, we can see that the FoM values obtained with 
the mock data are typically better than those for currently 
available data. This fact proves the capability of the forth- 
coming surveys to describe the evolutionary features of dark 
energy. 



5 CONCLUSIONS 

Parametrizing dark energy in a phenomenological way of- 
fers the possibility of progressing in the characterization of 
this main component of the Universe even though its origin 
has not been yet unveiled from a theoretical perspective. 
The extensive prior knowledge available on the topic hints 
that compliance with relevant astrophysical data favours 
parametrizations which are smooth and have two param- 
eters only. A second order requirement is restraint at high 
redshifts in the sense that dark energy should never red- 
shift slower than matter, but neither display an significant 
blueshift. 

Given those guidelines, in this work we present two new 
polynomial parametrizations of dark energy and we pay at- 
tention to some particular aspects which are important for 
model selection. On the one hand we examine the degree of 
correlation between the two parameters chosen in our pro- 
posal; specifically those are the EOS values at z = and 
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Figure 3. Total w and Icr error bands for our best fits for real (a) and mock data (b) compared with tiie results for the CPL 
parametrization (dashed lines and light contours). 



z — 0.5, which have been already shown to provide a conve- 
nient way to revisit and have already been shown to provide 
a preferred way to reconsider. On the other hand we resort 
to a genuinely Bayesian model selection criterion, the de- 
viance information criterion (DIG) for a better account of 
the improvements that our parametrizations represent. 

The astronomical tests we perform make use of SN lu- 
minosity data and the typical BAO related scales, so we are 
focusing on low redshift probes, although our parametriza- 
tions are suitable (as they are well behaved at high redshifts) 
for the analysis of CMB constraints which we hope to ad- 
dress in the future. 

The conclusions that emerge from our analysis is that 
our new parametrizations perform better in the sense that 
they allow to obtain tighter constraints on the dark energy 
EOS at present and its derivative, and they are also favoured 
by the statistical indicator we mention above (DIG). The 
main reason why we feel they fare better than their com- 
petitors is that they represent rather flexible perturbations 
of the ACDM scenario, which in many respects still remains 
the best description of the accelerated universe, and not sur- 
prisingly is often referred to as the concordance model. In 
contrast, either CPL or Wang are models which when taken 
as perturbations of AGDM are left with one free parameter 
only, and then one could expect less ability to accommodate 
themselves to the data. 

Note as well that the dataset we consider are not only 
some currently available, but we also simulate surveys to 
come, in particular, it is of interest our provision of syn- 
thetic measurements of the radial and transversal BAO 
scales as expected to be obtained with the EUCLID spec- 
troscopic survey, which as expected help us conclude that 
future surveys will decrease considerably our degree of igno- 



rance about dark energy evolution by providing far tighter 
constraints that may eventually lead to the conclusion that 
a cosmological constant is not the preferred candidate for 
the dark energy. 
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APPENDIX A: PARAMETER ESTIMATION 

The likelihood function, C{d\0,A4), is defined up to pro- 
portionality, as the probability of measuring the data d = 
{di, . . . ,d„} given the model M and its parameters taking 
the values = {^i, . . . ,9^} ( [Valiviita 2005 Denison et al. 
|2002| ). 

Despite our aim to keep the discussion in this section as 
general as possible, when we analyse particular datasets we 
will assume, as usual, that the measurements are normally 
distributed around their true value, so that 
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Figure 4. Confidence contours for the four parametrizations using current data. 
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-x^e)/2 
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The probability density function p{0\d,A4) of the pa- 
rameters to have values for the data, d, under the assump- 
tion that the true model is M is provided by Bayes' theorem 
( Valiviita|2005l ) 

C(d\e,M)n{0,M) 



p{6\d,M) 



J Cid\0,M)n{0,M)dO' 



(A2) 



where p{G\d,A4) and n{6,J\4) are the posterior and prior 
probability density functions (pdf) respectively (|Valiviita| 



2005[|Trotta|2004[|2007}|Cousins|1995}|Denison et al.|2002 | 

The prior pdf encodes all previous knowledge about the pa- 
rameters before the observational data have been collected. 
It can be regarded as a subjective procedure, but its use 
is compulsory in the Bayesian framework, which is the ap- 
proach used in theoretical frameworks where only one par- 
ticular realization of the measurement is available. 

Parameter estimation in the Bayesian framework is 
based on maximizing the posterior pdf p{6\d,A4), whereas 
in a "strict" frequentist approach one just maximizes 
jC{d\G,A4). When one uses fiat priors in the Bayesian ap- 
proach then the same conclusions are drawn from both ap- 
proaches and then the difference turns to be conceptual only 



(Trotta 2004 2007 Trottal. If the measured observables 



are independent form each other and Gaussian distributed 
around their true value, d{0), with a covariance matrix C, 



given by the experimental errors, maximizing C is equivalent 
to minimizing the chi-square function 



x\e) = d-'^ - did] 



d°'" - d(6») 



and for uncorrelated data Ca = Saa^ 



(A3) 



(A4) 



The second step toward constraining parameters sat- 
isfactorily is to construct credible intervals (Trotta 2004| 
which measure the degree of confidence that a certain data 
was generated by parameters belonging to the estimated in- 
terval, which gives the best fit parameters In the Bayesian 
approach, the credible intervals are drawn around the max- 
imum likelihood point, 6b f. After obtaining it by the mini- 
mization of the x^Wj the boundaries of the region contain- 
ing 100n% of likelihood are determined as the values of the 
parameters for which has increased by a certain quantity 



2 

with 



n = 1 



= 1 - 



2' 2 



(A5) 



(A6) 
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Figure 5. Confidence contours for the four parametrizations using simulated data. 



where T {i//2, A^^tfi) is the incomplete F function (Lazkoz 
et al.|2005| ), ( |Press et al.|1992| ). 

The la and 2a errors of the parameter 6i are given by 
the 68.30% and 95.45% credible interval contours, respec- 
tively. The upper limit is the maximum value of the contour 
and the lower one the minimum one. 



APPENDIX B: PIVOT COMPUTATION 
Bl Wang's parametrization 

The change of variables transforming the CPL parametriza- 
tion into Wang's is 



Wo + {1 — ac)Wa, 



(Bl) 



where Wc is the dark energy parameter EOS at the pivot 
redshift Zc for which some correlation related (good) feature 
is achieved. Wang's choice Wc = 0.5 leads to a low correla- 
tion situation, whereas we are letting Wc stand for the value 
associated with the lowest correlation possible, so our scope 
is more general. 

Using error propagation one gets 



(Jo 4- (1 — ac)aa + 2(1 — ac)ai2, 



(B2) 



directly from Eq. |B1[ In the remainder we will drop the 
absolute value given that we are only interested in a ac < 1 



situation. In contrast, if one solves Eq. |Bl| for Wa and then 
applies on it error propagation the result 



2 ctq + al — 2aoc 
a„ = 



(1 



follows. Combining our results above 
o"o + (1 ^ ac)a\2 = aoc 



(B3) 



(B4) 



is obtained, and then one deduces that total decorrelation 
[aoc = 0) is achieved for 



- (7o/o-12 



(B5) 



which corresponds for our results to Zc = 0.16,0.26 respec- 
tively for real and synthetic data. Note that according to our 
notation ai2 denotes in this case the non-diagonal element 
of the covariance matrix of wq and Wa, whereas aoc denotes 
the non-diagonal element of the covariance matrix between 
uiQ and Wc- 

But if one is just demanding a situation where aoc < o"i2 
it is very easy to deduce from Eq. |B4| that such condition is 
met whenever 



> ao/a\2. 



(B6) 
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B2 Chebyshev polynomial 

Provided luo.s ~ w\z=o.5, let us adopt for Wc the most gen- 
eral definition along the ideas discussed in their previous 
subsection. For our Chebyshev polynomial parametrization 
we have 

1 2 

Wc — — [—Gad^wo — 3wo.5 + 1) + ac{73wo — 63uio.5+ 

f 0) - 38wo + 45W0.5 - 4) . (B7) 

Following the same straightforward calculation as for the 
previous case we conclude there will be minimal correlation 
for 



^/729af2 - 2142(712^2 + IGSla^ + 63ai2 - 73crg 



12(3ai2 - 4cr2) 



(B8) 



which corresponds to Zc = 0.17,0.28 respectively for real 
and synthetic data. Sticking to our notation here ai2 denotes 
the non-diagonal element of the covariance matrix of wo and 
Wo. 5, and ctqc has the same meaning as in the case before. 

But again one can just want to meet the less restrictive 
requirement aoc < o"i2, which follows provided 

(3ac - 2)(6ac - 17)cri2 < (24a^ - rSa^ + 38) ctq . (B9) 



B3 Conventional polynomial 

Finally, following the same steps as before, but this time 
considering the conventional polynomial parametrization we 
have 



Wc = ;^(4(2 - 3ac){ac - 2)wo + 9{ac - 2){ac - l)iuo.5 + 



(2 — 3ac)ac + 3ac — 2). 



(BIO) 



The same route as for the two previous cases drives us to the 
conclusion that the minimal correlation situation is achieved 
for 



9(712 — 8(To 



(Bll) 



9(Ti2 - 12a2 ' 

which corresponds to Zc — 0.18, 0.29 respectively for real and 
synthetic data. Here the meaning of ai2 and aoc is exactly 
the same of the case before. 

Finally, one can deduce the condition for the less re- 
strictive requirement aoc < o"i2, to happen is simply 

^(7 - 3a,)(2 - 3a,)(7i2 < (2 - 3a,)(2 - ajcr^. (B12) 



APPENDIX C: EQUIVALENCE BETWEEN 

FOMdetf and FOMwang 



Following ( |Coe|2009| ) we get 



M^C"^M 



where 



(7i (Tl2 

2 

(Ti2 (T2 



(CI) 



(C2) 



and Mij = dpi/dp'j. In our case {pi,p2} ~ {wo,wo.5} and 
{Pi,p'2} = {wo,Wa}, and specifically wo.5 = mi + m2Wo + 



msWa with mi,m2,m3 coefficients different in the different 
parametrizations considered other than the CPL one, thus 



M : 



1 

m2 7713 



It is a simple matter of algebra to see that 



(Tl2 — m2(7i 



^ ' (Tl2 — m2(Ji m%ol — 2m2(Tl2 + (^2 



(C3) 



(C4) 



ma 



7713 



and \/ det C = Vdet C /m^. Therefore we get the simple 
relationship between the two definitions of FoM 

FoMdETF = 77l3FoMwang. (C5) 
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